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Abstract 

Background: Frog virus 3 (FV3) is the type species of the genus Ranavirus, and in the past few decades, FV3 
infections have resulted in considerable morbidity and mortality in a range of wild and cultivated amphibian 
species in the Americas, Europe, and Asia. The reasons for the pathogenicity of FV3 are not well understood. 

Findings: We investigated three FV3 isolates designated SSME, wt-FV3, and aza-C r , and reported that our wt-FV3 
and aza-C strains showed similar levels of virulence, while SSME was the least virulent in an in vivo study with 
Lithiobates pipiens tadpoles. Using 454 GS-FLX sequencing technology, we sequenced SSME and compared it to the 
published wt-FV3 genome. SSME had multiple amino acid deletions in ORFs 49/50L, 65L, 66L, and 87L, which may 
explain its reduced virulence. We also investigated repeat regions and found that repeat copy number differed 
between isolates, with only one group of 3 isolates and 1 pair of isolates being identical at all 3 locations. 

Conclusions: In this study we have shown that genetic variability is present between closely related FV3 isolates, 
both in terms of deletions/insertions, and even more so at select repeat locations. These genomic areas with 
deletions/insertions may represent regions that affect virulence, and therefore require investigation. Furthermore, 
we have identified repeat regions that may prove useful in future phylogeographical tracking and identification of 
ranaviral strains across different environmental regions. 
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Background 

Global amphibian populations have declined consider- 
ably in recent years, in part due to habitat fragmenta- 
tion, pollution, and the chytrid fungus Batrachochytrium 
dendrobatidis [1,2]. More recently, certain members of 
the family Iridoviridae have also been associated with 
amphibian decline. The family Iridoviridae is comprised 
of large, cytoplasmic, double stranded DNA viruses with 
icosahedral capsids [3], and is divided into five genera: 
Iridovirus, Chloriridovirus, Lymphocystivirus, Megalocyti- 
virus, and Ranavirus [4]. Specifically linked to amphib- 
ian mortality in this family are infectious diseases caused 
by members of the genus Ranavirus. In past research, 
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ranaviruses received little attention as most infections 
were deemed subclinical; however, recent ranavirus in- 
fections have resulted in considerable morbidity and 
mortality in a range of wild and cultivated amphibian 
species in the Americas, Europe, and Asia [5-8]. It has 
been reported that 43% of known amphibian die-offs in 
the USA from 2000 to 2005 were due to ranaviruses [9], 
and that from 1996-2001 ranaviruses were isolated from 
most of the amphibian mortality events in North America 
[10]. Detection of these outbreaks could be due to better 
surveillance, increased environmental awareness, the mu- 
tation of viral species creating highly pathogenic strains, 
or environmental changes resulting in host immune sup- 
pression [11], Ranaviruses have become a significant cause 
of death and disease in amphibians, and thus investigation 
into these viruses is warranted from a virological, com- 
mercial, and ecological standpoint [11]. 
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Frog virus 3 (FV3) is the type species of the genus 
Ranavirus [4]. FV3's genome is 105,903 base pairs (bp) 
comprised of 98 open reading frames (ORFs) [12]. De- 
pending on factors such as strain virulence and host im- 
mune response, infection with FV3 may or may not lead 
to mortality. However, in susceptible amphibians, FV3's 
necrotic and apoptotic effects cause systemic, chronic 
cell death in multiple internal organs, resulting in death 
of the host within a few days to several weeks [2,13,14]. 
FV3 infection is also marked by cutaneous signs, includ- 
ing ulceration of the skin, and erythema and swelling of 
the limbs and body. In fatal cases, intracoelomic lesions 
are often present, including haemorrhages of the kidneys 
and reproductive organs, and pale, swollen livers [15]. 
While our understanding of ranavirus pathogenicity has 
improved over the last decade, there is still a need for 
the research community to more fully describe the de- 
terminants of virulence variation. Elucidation of this area 
will likely heavily rely on the genetic analysis and com- 
parison of ranaviruses that differ in host range and 
virulence. 

Genetic comparison of related DNA viruses has proven 
to be an important tool in classifying viral strains and 
understanding the epidemiology and evolution of differ- 
ent genotypes [16]. By analyzing genetic differences, re- 
searchers can link clinically significant alterations with 
molecular changes, and better understand viral origins 
and evolution. Increasingly, new viral strains are being 
identified based on the systematic analysis of sequence 
data, including short amino acid insertions, transla- 
tional stop codons, and single amino acid deletions 
[17]. Analysis of viral isolates has led to the discovery 
of new viral genotypes, as well as better understanding 
of the functional genetic differences among strains [17]. 
For instance, the entire genome of a virulent strain of 
duck entiritis virus was recently sequenced and com- 
pared to the genomes of an attenuated strain and an- 
other virulent strain [18]. The results indicated several 
nucleotide insertions/deletions and frame-shift mutations 
effecting ORF initiation or termination [18]. These find- 
ings allowed the researchers to identify possible virulence 
factors and provided information on ORFs that are chan- 
ged during serial passage. 

In addition to nucleotide insertions/deletions, variation 
between viral genomes may occur at repeat regions. 
Eaton et al. [19] identified repetitive sequences in the ge- 
nomes of various ranaviruses with high copy number 
variation. Repetitive sequences are commonly classified 
into one of three groups: macro, mini, or micro satellites 
[20]. The repeats that we will examine in this study con- 
tain less than 400 bp comprised of 9-19 bp repeating 
units. Thus, we suggest that these repeats are some vari- 
ation between micro and mini satellites, and for the pur- 
poses of our analysis, will be referred to as short tandem 



repeats (STRs). A selection of these repeat regions will 
be used to analyze FV3 isolates in order to further inves- 
tigate the fine scale, genetic differences present in vari- 
able regions. 

Past studies on the genetic variation between DNA 
viral strains have allowed for the detection of minute 
genetic changes that would otherwise have gone un- 
noticed. These kinds of changes have proved useful in 
explaining phenotypic differences and evolutionary his- 
tories. The purpose of the present investigation is to nar- 
row our focus even further by comparing the genomes 
of closely related FV3 isolates, including those with vary- 
ing levels of virulence. This will be done in an attempt 
to explore the genetic diversity present in strains of FV3, 
with the ultimate goal of further elucidating the possible 
genetic basis behind FV3's unpredictable infectious 
behaviour. 



Results 

FV3 strains differ in virulence during in vivo infection 

In order to determine if FV3 strains SSME, wt-FV3, and 
aza-C r induced different degrees of infection, L. pipiens 
tadpoles were exposed to each strain. Tadpoles were 
exposed to FV3-infected water for 12 h before being 
transferred to FV3-free water for the remainder of the ex- 
periment. The tadpoles were monitored for 41 days, at 
which time all tadpoles had either died or reached meta- 
morphosis. The results of the survival analysis revealed 
significant differences in strain effect on tadpole mortality 
rates (X 2 = 21.3, p < 0.01). Although deaths triggered by 
wt-FV3, aza-C r , and SSME stocks peaked between 15 and 
18 days post infection, overall mortality was higher follow- 
ing infection with wt-FV3 and aza-C r virus (>96%) than 
with SSME (84%). Wt-FV3 induced a mortality of 97%, 
and aza-C r a mortality rate of 96% (Figure 1). In contrast, 
SSME infected tadpoles reached a mortality plateau after 
28 days, and showed a mortality rate of 84%. The control 
tadpoles had the lowest mortality, with a rate of 66% after 
41 days (Figure 1). These data demonstrate that the 
closely related FV3 strains, wt-FV3 and aza-C r [21,22], are 
associated with higher mortality than SSME (Figure 1). Of 
concern was the high mortality seen even in control tad- 
poles (Figure 1). We did confirm that the control tadpoles, 
both those that survived and those that died, were rana- 
virus negative (data not shown). This mortality could be 
the result of overcrowding of tadpoles during infection. In 
addition, this mortality is occurring during the first few 
days of the experiment at a time when the tadpoles are 
very small and very weak, very sensitive to stressful condi- 
tions. At the beginning of the experiment, the tadpoles are 
transferred from tank to tank that may account for the 
high mortality rate early on in the experiment or may have 
resulted in high mortality of a sublethal viral dose. 
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Figure 1 Variable virulence of three FV3 strains during tadpole infection. Survival analysis of L pipiens tadpoles over 41 days of infection 
with aza-C r (square), SSME (circle), wt-FV3 (triangle), and Control (diamond). Tadpoles were exposed to FV3 strains for 12 hours, and then along 
with the virus-contaminated water were transferred to a tank of dechlorinated water for the remainder of the experiment. Survival analysis and 
failure time analysis was done following the Kaplan & Meier product limit method associated with Chi square and Gehan's Wilcoxon tests [37]. 



Genomic sequencing of SSME 

In order to better understand the possible genetic basis 
of phenotypic variation among strains, we sequenced the 
SSME strain using 454 GS-FLX technology. Our se- 
quenced genome showed high similarity to the published 
FV3 sequence (wt-FV3) [12], with a sequence identity of 
9879% and average genome sequencing coverage of 51x. 
Results revealed that SSME differed from the wt-FV3 
genome length of 105,903 bps and ORF number of 98; 
instead, SSME had a length of 105,070 bps, and a total 
of 95 predicted ORFs (Table 1). Despite high overall se- 
quence identity between SSME and FV3, marked differ- 
ences were noted in several regions of the genome 
(Table 2). 

Gene variation between ranaviruses 

In order to investigate the possible genetic variation 
within coding regions and regions of high variability be- 
tween closely related viral isolates, we compared our se- 
quenced SSME genome with four published ranavirus 



genomes. The genomes compared included: SSME, wt- 
FV3 [17], rana grylio iridovirus (RGV) [23], soft-shelled 
turtle iridovirus (SSTIV) [24], and tiger frog virus (TFV) 
[25]. We began by focusing our analysis on predicted 
coding regions, as variation in these areas may have con- 
sequences in pathogenicity or the kind of disease a virus 
causes. 

Upon analysis of selected ranaviral genomes, a 757 bp 
deletion was discovered in SSME, deleting the entire 65L 
coding region and the majority of the 66L coding region 
(Figure 2). This deletion was not present in wt-FV3 or in 
related ranaviruses RGV, SSTIV, and TFV. However, se- 
quence alignments indicated reorganization of 66L in 
RGV, SSTIV, and TFV, with a 139 bp insertion found in 
their corresponding 66L regions (Figure 2). Therefore, 
there are at least 3 different genomic presentations of this 
region in ranaviral genomes, making it highly variable 
among isolates. 

Another coding region demonstrating variation was 
50L. An in frame deletion of 13 amino acids was found 
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Table 1 Description of nucleotide start/stop locations and 
amino acid length of ORFs in SSME, along with 
predicted functions 



Table 1 Description of nucleotide start/stop locations and 
amino acid length of ORFs in SSME, along with 
predicted functions (Continued) 



ORF 



Start/stop (#aa) 
SSME 



Predicted function/conserved domains 



1R 272-1042 

255aa 
2L 1649-2629 

325aa 
3R 3436-4752 

437aa 
4R 4793-4975 

59aa 

5R 5408-6022 

203aa 
6R 6025-6252 

74aa 

7R 7043-7429 
127aa 

8R 7521-11402 

1292aa 
9L 11771-14617 

947aa 

10R 14633-15046 
136aa 

11R 15396-15608 
69aa 

12L 15674-16567 
296aa 

13R 17108-17314 
67aa 

14R 17329-17688 
118aa 

15R 17784-18752 
321aa 

16R 19032-19859 
274aa 

17L 20100-21608 
501aa 

18 L 21645-21881 TM 
77aa 

19R 21933-24515 
859aa 

20R 24562-25008 
147aa 

21 L 25231-25890 
218aa 



Replication factor &/or DNA packing protein 



Myristylated membrane protein, 
DUF230 poxvirus protein, TM* 

IIV6 ORF 229L, SP 



TM 



FPV ORF 250, US22, herpes virus early 
nuclear protein 

SP* 



SP 



DNA-dependent RNA polymerase largest 
subunit 

NTPase, SNF2 family, N-terminal, 
helicases C-terminal, DEAD/H helicases 

TM 



TM 



Unknown protein, SP 



SP 



Unknown protein, SP 



AAA-ATPase, poxvirus A32 



Integrase homologue 



SP 



Similar to LCDV1 or f10L conserved 
uncharacterized protein, SP 

Unknown protein, TM 



ISKNV ORF 56L-like protein 



22R 



23R 



24R° 



25R 



26R 



27R 



28R 



29L 



30R 



31R 



32R 



33R 



34R 



35L 



36L 



37R 



38R 



39R 



40R 



41 R 



42L 



26020-28941 
972aa 

29319-30467 
381 aa 

30877-31947 
355aa 

32141-32929 
261 aa 

32996-33226 
75aa 

33757-36669 
969aa 

36718-37206 
161aa 

37385-37681 
97aa 

37883-38035 
49aa 

38097-38516 
138aa 

38566-40458 
629aa 

40541-40732 
62aa 

40875-41195 
105aa 

41287-41748 
152aa 

41761-42384 
206aa 

42778-43413 
210aa 

43554-45251 
564aa 

45357-45707 
115aa 

45795-46343 
181aa 

46725-50222 
1164aa 
50718-50975 
84aa 



D5 family NTPase, ATPase 



SP 



SP 



P31K protein 



Truncated elf-2a homologue 



Tyrosine kinase, CAP 10, putative 
lypopolysaccharide modifying enzyme 

Unknown protein, SP 



SP 



TM 



SP 



Neurofilament triplet H1 protein 



Unknown protein, TM 



Human parainfluenza virus 1 L like protein, TM 



SP 



SP 



NIF/NLI interacting factor 



Ribonucleoside diphosphate reductase 
alpha subunit barrel domain 

Hydrolase of the metallo-beta-lactamase 
superfamily 

TM 



RRV ORF-2-like protein, SP 



SP 
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Table 1 Description of nucleotide start/stop locations and 
amino acid length of ORFs in SSME, along with 
predicted functions (Continued) 



43 R 


50974-51491 




1 71 aa 


44 R 


51512-51697 




60aa 


45 L 


51973-52383 




1 35aa 


46L§ 


52437-53000 




186aa 


47L 


53125-53541 




137aa 


48 L 


53544-53795 




82aa 


49/50L1 


53904-55451 




514aa 


51R 


55531-57216 




560aa 


52L 


57473-58540 




354aa 


53R 


58878-60446 




521aa 


54L 


60677-60907 




75aa 


55L 


60945-62240 




430aa 


55R 


61089-62228 




378aa 


56R 


62328-62765 




144aa 


57R 


62879-64375 




497aa 


58R 


64700-65413 




236aa 


59L 


65964-67022 




351aa 


60R 


67184-70225 




1012aa 


61L 


70234-70416 




59aa 


62L 


70859-74524 




1220aa 


63R 


74903-75397 




163aa 



Table 1 Description of nucleotide start/stop locations and 
amino acid length of ORFs in SSME, along with 
predicted functions (Continued) 



TM 64R 

SP 65 L # 

66L # 

LCDV1 ORF-88-like protein, SP 



RGI 47L-like protein, SSTIV 049L-like 
protein, Neurofilament triplet HI -like 
protein, microneme/rhoptry antigen 

SP 



SP 



3-beta-hydroxy-delta 5-C27 steroid 
oxidoreductase-like protein, TM 

LCDV1 ORF-20-like protein, SP 
Nuclear calmodulin-binding protein 
Helicase-like protein 
FV3 40-kDa protein, SP 
SP 



SP 



RGV 9807 unknown protein, SP 



SP 



DNA-directed RNA polymerase II second 
largest subunit RNA polymerase domain 
6, 7, 3, 2 beta subunit 



67L 



68R 



69R 



70R 



LCDV1 ORF 58-like protein, RGI ORF 

50L-like protein, SSTIV ORF 052-like 71 R 

protein, SAP DNA binding domain 

Unknown, SP 



72L 



73L 



74L 



75L 



76R 



77L 



78L 



Phosphotransferase, S-TKc, Serine/ 

Threonine protein kinase 79R 



80L 



81R 



DNA polymerase-like protein, DNA 

polymerase family B exonuclease 82R 



83R 



dUTPase-like protein 



85 R 



75522-75809 

94aa 

N/A 

76106-76157 
16aa 

76212-77375 
386aa 

77657-77944 
94aa 

78080-78346 
87aa 

78364-78738 
123aa 

78778-7901 1 
76aa 

79068-79784 
238aa 

80232-81206 
323aa 

81381-82493 
369aa 

82525-82779 
83aa 

82842-83063 
72aa 

83060-83407 
114aa 

84007-84645 
21 laa 

84781-86499 
571aa 

87122-88237 
370aa 

88293-88571 
91aa 

88700-89173 
156aa 

89623-90267 
213aa 

90652-91389 
244aa 

91464-92051 
194aa 



lnterleukin-1 beta convertase precursor, 
Caspase-recruitment domain/DEATH 

SP 
SP 



Ribonucleoside-reductase diphosphate 
beta subunit-like protein 

SP 



Unknown protein, TM 



SP 



SP 



SP 



NTPase/helicase-like protein 



SP 



LITAF/PIG7 possible membrane associated 
motif in LPS-induced tumor necrosis factor 
alpha factor, TM 

SP 



LCDVORF 2-like protein, SP 



SP 



ATPase-dependent protease 



Ribonuclease Ill-like protein 



Transcription elongation Factor Sll, 
C2C2 zinc finger 

Immediate-early protein ICP-18 



Cytosine DNA methyl-transferase 



LCDV 1-like proliferating cell nuclear 
antigen, SP 

Deoxynucleoside kinase 
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Table 1 Description of nucleotide start/stop locations and 
amino acid length of ORFs in SSME, along with 
predicted functions (Continued) 



Table 2 Genetic variation present within the SSME 
genome as compared to the wt-FV3 reference genome 



86L 



87L* 



89R 



90R 



91R 



92R 



93L 



94L 



95R 



96R 



97R 



98R 



92441-92626 
60aa 

92979-94778 
598aa 

94811-95263 
149aa 

95331-96497 
388aa 

96590-97981 
462aa 

98105-99292 
394aa 

99627-99872 
80aa 

100054-100221 
54aa 

100331-100798 
154aa 

100891-101982 
362aa 

102784-103455 
222aa 

103538-103951 
136aa 

104649-104849 
65aa 



SP 



Unknown protein, SP 



Evrl-air-augmenter of liver regeneration 



SP 



Major capsid protein 



Immediate-early protein ICP-46 



SP 



SP 



Regina ranavirus P8.141 C-like protein, TM 



DNA repair protein RAD2, Xeroderma 
pigmentosum G N-region, Helix-hairpin-helix, 
Class 2 (Pol I) family 

Unknown protein, SP 



Myeloid cell leukemia protein, MCL-1,TM 



SP 



*Note: SP-surface protein, TM-transmembrane domain. 

°Loss of original, predicted start codon in SSME; new start codon 

27 bp downstream. 

§Loss of original, predicted stop codon in SSME; new stop codon 
319 bp upstream. 

1139 bp deletion in original 50L ORF of SSME. Loss of ORF 50L's original, 
predicted stop codon in SSME; ORFs 49L and 50L are now combined as 
49/50L. 

#Deletion of 65L ORF and majority of 66L ORFs in SSME. 
**18 bp deletion in 87L ORF of SSME. 

in 50L when compared to wt-FV3, present only in the 
SSME genome (Figure 3). In addition, 50L had a single 
nucleotide deletion resulting in a frameshift mutation 
in SSME, RGV, and SSTIV when compared to wt-FV3. 
This frameshift mutation resulted in the loss of 50Ls 
stop codon, which led to the combination of 49L and 
50L into a single ORF. These two ORFs merged in 
frame, creating 49/50L. In contrast, TFV did not have 
this single nucleotide deletion, nor did it display the 
merger of its 49L and 50L equivalent ORFs (Figure 3). 
This is consistent with TFV being the most divergent 
virus of the group. 



Region 

21,933-24,515 bp 
30,851 bp 



51,102 bp 
52,769 bp 



54,805 bp 

54,938 bp 
76,105 bp 

94,712 bp 
10,4050 bp 



Change 



Affect 



24 bp substitutions 24aa substitutions & 9aa 

^ . . . insertion in 19R 

27 bp insertion 



1 bp substitution 



1 bp deletion 
1 bp deletion 

1 bp deletion 

39 bp deletion 
757 bp deletion 

18 bp deletion 
67 bp deletion 



Loss of a start codon in 24R; 
new start codon is 27 bp 
downstream, thus 
shortening ORF by 9aa 

Frameshift mutation in 43R 

Frameshift mutation and loss of 
original stop codon in 46L; 
new stop codon is 319 bp 
upstream, thus lengthening 46L 

Loss of original stop codon in 
50L, leading to the combination 
of ORFs 49L and 50L 

13aa deletion in 49/50L 

Deletion of all of 65L and the 
majority of 66 L 

6aa deletion in 87L 

67 bp deletion in non-coding 
region between 97R and 98R 



Genome region involved, nucleotide change, and end affect of the variation 
is presented. 



Further ranaviral gene variation was discovered in 43R 
in the form of a single nucleotide deletion, resulting in a 
frameshift mutation. Similar deletions were found in 
SSME, RGV, SSTIV, and TFV, leading to pronounced 
differences in amino acid sequences when compared to 
wt-FV3 (Figure 4). It is worth noting that although 
SSME, RGV, and SSTIV share the same frameshift muta- 
tion the genes are not identical as there are several sin- 
gle amino acid changes that lead to modest variation 
(Figure 4). In contrast, despite having a similar nucleo- 
tide deletion, TFVs proposed amino acid sequence was 
markedly different from other analyzed ranaviruses. An- 
other single nucleotide deletion was found in 46L 
(Tables 1 and 2), which led to the loss of the original 
stop codon and the extension of the ORF by 319 bps 
(Tables 1 and 2). This was found in all genomes other 
than TFV and wt-FV3. These data demonstrate that 
even between very similar viral genomes, there is consid- 
erable gene variation, and that when viruses differ in 
pathogenicity, it could be related to variation in coding 
regions. 

Repeat regions can vary in copy number between 
viral isolates 

Genomic sequencing, and subsequent analysis, of SSME 
in comparison to wt-FV3, provided a general overview 
of the genetic variation between closely related FV3 iso- 
lates, particularly in coding regions. We then decided to 
explore possible genetic variation in highly variable 
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wt-FV3 



64R 





65L 


66L 






SSME 



RGV, SSTIV & 
TFV 









64R 









64R 




65L 






66L 













Figure 2 65L & 66L deletion in SSME. Alignment of the 64R-66L region across ranaviruses. A 757 bp deletion is present in the SSME genome, 
spanning from 76,1 13-76,869 bp in wt-FV3, deleting all of 65L and most of 66L A 139 bp insertion is also present in the 66L homologous 
regions of RGV, SSTIV, and TFV genomes. 



repeat regions, and thus define these genetic regions with 
greater inter-strain variability. To achieve this, we chose 
three repeat regions (suspected to be polymorphic in copy 
number) to investigate. These regions included: Region 1 
(22,499-22,574 bps); Region 2 (52,443-52,747 bps); and 
Region 3 (54,948-54,986 bps), based on the wt-FV3 se- 
quence. We performed a repeat analysis on our sequenced 
SSME genome that we had previously evaluated for cod- 
ing region variability, and reference genomes RGV, SSTIV, 
and TFV. We also decided to analyze the repeat regions of 
wt-FV3 and aza-C r in order to check for repeat number 
stability between the two strains over the multiple viral 
passages they have undergone since the creation of aza-C r 
in 1987 [22]. In addition, we expanded our analysis to in- 
clude 6 environmental samples isolated from the same 
waterway in Manitoulin, Kagawong, ON, Canada. These 



wt-FV3 



RGV & SSTIV 



SSME 



TFV 









50L 







49/50L 



49/50L 



Al3aa 



Figure 3 Variation in ORF 49/50L. Sequence alignment of the 
49/50L region across ranaviruses. A 39 bp deletion is present in the 
SSME genome, spanning from 54,909-54,947 bp in wt-FV3. A single 
nucleotide deletion is present in SSME, RGV, and SSTIV, leading to 
the loss of 50L's original stop codon and the merger of 49L and 50L 
into one ORF. Different colours are used to represent the portions of 
TFV's ORFs that are homologous to parts of 49L and 50L respectively. 



isolates were designated as: E3, E4, E5, F4, F6, and G4; a 
lab FV3 sample (ATCC VR-567™) was also used [26]. 

The information acquired through analyzing these re- 
peat regions indicated that variation in repeat copy num- 
ber is present between ranaviral isolates, including those 
from one geographic location (Figure 5). In the first re- 
gion, repeat copy number ranged from 2-5 copies, with 
2 copies being the most common and 5 copies being 
unique to isolate F6. Wt-FV3 and aza-C r shared a dis- 
tinctive repeat sequence, containing a G rather than an 
A at the 9 th nucleotide position (Figure 5). Region 2 
proved to be more variable than region one, with repeat 
copy number ranging from 4-18 copies, while Region 3 
was found to consist of two repeat sequences. These 
Region 3 repeats varied in copy number between iso- 
lates, with a maximum copy number of 4. However, in 
some instances, the first sequence of nucleotides charac- 
teristic of this region did not repeat (Figure 5). 

Although copy number variation was present amongst 
isolates, when samples were compared across all three 
repeat regions, some were found to be identical; such 
isolates included F4, E5, and E3 (Figure 5), which were 
isolated from the same geographic region. In addition, 
isolates wt-FV3 and aza-C r were also found to be identi- 
cal across repeat regions (Figure 5), which is consistent 
with their common origin [22], This data taken as a 
whole suggests that viral isolates, even those from one 
geographic location, can display both variation and iden- 
tical patterns across repetitive regions of the genome. 
Thus, these repeats could act as a fingerprint to discrim- 
inate between viral isolates. 

Discussion 

Our analysis of related ranaviruses presents a novel ap- 
proach to genomic comparison that differs from other 
studies. We analyzed the genomes of closely related iso- 
lates of FV3 through 454 GS-FLX technology and STR 
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SSME MEEVDGCAGPNS EAGALTAGALTAGAFAVTAGAG VAGAGVAGVAGVAGAA — GAA 

RGV MEEVDGCAGPNS EAGALTAGALTAGAGAAAAGAGAGV- AVAGVAVAGV AVAGVA 

SSTIV MEEVDGCAGPNS EAGALTAGALTAGAGAAVAGAVAGV- AV AGV AVAGVA 

TFV MEEVDGCAGPNSEAGALTAGSLTAVALTAGVDGACPKDGKTEVGLKGVNPVKAGAFAVTA 

Wt-FV3 MEEVDGCAGPNS EAGALTAGALTAGAFAVTAGAGVAGAGVAGVGWCSWCSWCSWCWCSWC 



bbMt GAGVAGAAGAAGAAGAGAAGAGAAGAGAAGAGAA GAGA- AGAGAAGAC 

RGV VAGVAVAVAGAGVAGVAVAGAGVAGVAVAGAGAAGAGVAVAGVAVAGAGAAGAAGVAGAC 

SSTIV VAGVAVAVAGAGVAGVAVAGAGVAGVAVAGAGAAGAGVAVAGVAVAGAGAAGAAGVAGAC 

TFV GAGV DGAGVDGAGVDGAGVDGAGVSGAGVS GAGV- S GAGVSGAC 

Wt-FV3 SWCWC-SWCWCSW CW CS WCWCSWCWCSVJCW-CSWCWCSWC 



55ivit P-KDGKTEVGLKGV74PVKAGAFAVTAGAGVAG--AGVAGAAGACPKDGKTEVGLKGVNPV 

RGV P-KDGKTEVGLKGVNPVKAGAFAVTAGAF AVTAGAGVAGACPKNGKTEVGLKGVNPV 

SSTIV P- KDAKT EVGLT GVN P VRLVP LQS LLVP LQS LLVLV* LVP VQ RMP DRG GLDGCKSC 

TFV P-KDGKTEVGLKGVNPVKAGAFAVTA CAGACPKDGKTEVGLKGVNPV 

Wt-FV3 LSKGWEDRGGLEGCKSCKGWCLCSHCWCWCSWC — WCSWCSWCLSKGVIEDRGGLEGCKSC 



SSME KAGAFAVTAGAGALT 

RGV KAGAGALTAGAGALT- - - 

SSTIV k GWC * AS LLVL * Q * L 

TFV KAGA 

Wt-FV3 K GWCLCSHCRCWSIN 

Figure 4 Variation in 43R ORF. Amino acid sequence alignment of the 43R ORF in ranaviruses. A single nucleotide deletion in SSME, RGV, SSTIV, 
and TFV led to frameshift mutations, causing significant amino acid variability amongst genomes. 



comparison. The scale (fully sequenced genomes), and 
the nature of comparison (using viral isolates of FV3), 
set our investigation apart from past studies. We found 
that the 3 strains we examined displayed slightly differ- 
ent levels of virulence during in vivo studies. By sequen- 
cing the one genome yet to be sequenced, we were able 
to highlight areas that may be important in generating 
infectious phenotypes of FV3. This kind of analysis has 
never been done for FV3 or related ranaviruses: thus, it 
provides greater insight into the genetic variation among 



these closely related DNA viruses and the possible gen- 
etic basis of ranaviral virulence. 

Our results demonstrate that genetic variation is present 
between closely related FV3 isolates in both coding and 
non-coding regions. The SSME genome was sequenced 
and compared to the published wt-FV3 genome, along 
with related ranaviral genomes RGV, SSTIV, and TFV. 
Comparisons revealed that SSME was divergent from wt- 
FV3, and aza-C r . This variation could be due to the fact 
that SSME was isolated from a spotted salamander rather 



Viral Isolate 


Region 1 


Region 2 


Region 3 


F4, E5, &E3 


x2 


xl7 


x2 


x4 


- CTGCGAGAAAGTCTC - 


- CTTACGCTTCACAGGGGA- 


- TCAGGCTGC 


TCAGGCTGT- 


E4 


x2 


xl5 


x2 


x4 


- CTGCGAGAAAGTCTC - 


- CTTACGCTTCACAGGGGA- 


-TCAGGCTGC 


TCAGGCTGT- 


G4 


x4 


xl4 


x2 


x4 


- CTGCGAGAAAGTCTC - 


- C TTACGC TT C AC AGGGGA - 


-TCAGGCTGC 


TCAGGCTGT- 


wt-FV3 and aca-C r 


x3 


xl6 




x2 


- ctgcgagaa|gtctc - 


- CTTACGCTTCACAGGGGA- 


-TCAGGCTGC 


TCAGGCTGT- 


SSME 


x3 


xl6 




x2 


- CTGCGAGAAAGTCTC - 


- C TTACGC TT C AC AGGGGA - 


-TCAGGCTGC 


TCAGGCTGT- 


F6 


x5 


xl8 


x2 


x2 


- CTGCGAGAAAGTCTC - 


- CTTACGCTTCACAGGGGA- 


-TCAGGCTGC 


TCAGGCTGT- 


VR-567 


x4 


x4 


x2 


x2 


- CTGCGAGAAAGTCTC - 


- C TTACGC TT C AC AGGGGA - 


-TCAGGCTGC 


TCAGGCTGT- 


RGV 


x2 


xlO 






- CTGCGAGAAAGTCTC - 


- CTTACGCTTCACAGGGGA- 


-TCAGGCTGC 


TCAGGCTGT- 


SSTIV 


x2 


xl8 






- CTGCGAGAAAGTCTC - 


- C TTACGC TT C AC AGGGGA - 


-TCAGGCTGC 


TCAGGCTGT- 


TFV 


x3 

- CTGCGAGAAAGTCTC - 


- TGAAACAGACTTGCTGTCTC - 


- TCAGGCT - CCACAACAGCCT - 



Figure 5 STR variation summary. Short tandem repeats (STRs) vary between viral isolates at three tested locations. Results reveal three isolates 
(F4, E5, and E3), and two isolates (wt-FV3 and aza-C r ), that have the same STR copy number at each location. Areas highlighted in blue are full 
repeats, with their coinciding copy number above. Areas highlighted in green represent a sequence that may/may not repeat but is very similar 
to the STR that follows it in the genome, and areas highlighted in pink represent nucleotides that differ from the STR found in the majority of 
ranaviruses analyzed. 
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than anurans, and so the strain may have evolved in 
order for it to better adapt to its novel host, as is seen 
during serial passage [27]. For example, the pathogen- 
icity of Dengue virus was altered by serial passaging 
Dengue virus 27 times which resulted in 25 nucleotide 
changes between 2 strains [27]. This finding is interest- 
ing given the fact that our results also showed that 
SSME was the least virulent of our 3 strains during an 
in vivo study in anurans (Figure 1). Thus, genetic muta- 
tions observed in SSME, such as deletions in ORFs 65L, 
66L and 49/50L, could have had an effect on strain 
virulence in anuran hosts. Also, past sequencing results 
[12] as well as our own genomic sequencing of wt-FV3 
and aza-C r genomes revealed only 13 nucleotide differ- 
ences (data not shown) between the strains, supporting 
our in vivo finding that the two strains cause essentially 
the same amount of tadpole mortality (97% and 96% 
respectively). 

One of the most significant genetic variations found 
in SSME was a 757 bp deletion that deleted all of the 
65L coding region and most of the 66L coding region 
(Figure 2). Among other ranaviruses, the functions of 
these genes have not been determined. Further sequen- 
cing of ranaviruses has shown that 65L is present in 
both RGV and SSTIV as 69L and 68L respectively 
[23,28], while this form of 66L appears to be unique to 
FV3, as other ranaviruses have a 139 nucleotide inser- 
tion in their related 66L regions (Figure 2). Another 
genetic variation specific to SSME was a 13 amino acid 
deletion in ORF 50L (Figure 3). Interestingly, a second 
variation was found in this area in the form of single 
nucleotide deletion. This deletion lead to the merger of 
ORFs 49L and 50L in all examined ranaviruses other 
than TFV and wt-FV3 (Figure 3). In terms of function, 
49L has multiple SAP motifs, which are DNA/RNA 
binding domains predicted to be involved in chromo- 
somal organization and DNA replication [29]. Thus, 
this new 49/50L ORF may function in viral replication. 

Single nucleotide deletions were found in multiple rana- 
viruses within the 43R genes (Figure 4) and 46L (Tables 1 
and 2) genes. In 43R, the deletion was present in all rana- 
viruses analyzed other than wt-FV3, and resulted in a 
frameshift mutation (Figure 4). In 46L, the original stop 
codon was lost, leading to the extension of the ORF by 
319 bps (Tables 1 and 2). This was found in all genomes 
other than TFV and wt-FV3. Supposedly, 46L encodes for 
a neurofilament triplet HI -like protein [12]. However, the 
extended version of 46L that we discovered has a putative 
conserved domain known as a microneme/rhoptry anti- 
gen in the area previously thought to be non-coding. 
Micronemes and rhoptries are organelles possessed by 
Apicomplexa protozoans that secrete proteins involved in 
parasite entry into a host cell, specifically possessing 
protein-binding motifs that recognize ligands on the host 



cell surface [30]. Although usually associated with proto- 
zoan parasites, these microneme/rhoptry antigens found 
in 46L could give further indication as to 46Ls function. 

The changes present in 43R, 46L, 49/50L, 65L, and 
66L represent the main variable regions amongst the re- 
lated ranaviral isolates we analyzed. This suggests that 
these are changeable areas across ranaviruses, and could 
be used in the future to help explain variable infectious 
phenotype. Moreover, multiple amino acid deletions 
present in 65L, 66L, and 49/50L, were limited to the 
SSME strain, which displayed the lowest level of viru- 
lence during tadpole infection. Thus, changes that are 
unique to SSME may present areas of the genome that 
are particularly effective in viral attenuation, specifically 
in an alternative host. As FV3 mortality and morbidity 
continues to worsen and fluctuate across environmental 
regions, examination of these genomic areas may prove 
useful as an initial way to investigate the genetic basis 
behind infective changes [6]. Further research could be 
used to explain variations in the virulence of different 
FV3-like isolates. 

As we had already identified variation in the coding re- 
gions of closely related FV3 isolates, we decided to further 
our understanding of variation within highly variable sites 
by investigating 3 recently identified STR regions [19]. 
These are known to be variable areas: for instance, al- 
though FV3 and SSTIV share 99% genome sequence iden- 
tity, they share only approximately 50% of repeats in 
common [19]. Thus, we predicted that these repeat re- 
gions would have greater inter-strain variability that would 
provide useful information when trying to understand 
overall genetic variability between ranaviral isolates. In 
order to test this prediction, we analyzed our sequenced 
SSME genome and reference ranaviral genomes, along 
with 6 environmental samples isolated from the same 
waterway. We also sequenced wt-FV3 and aza-C r to check 
for repeat number stability across their past viral passages. 
Analysis revealed that repeat copy number was variable 
between isolates, even between those from the same 
geographic location, but that there was some conserva- 
tion (Figure 5). Specifically, wt-FV3 and aza-C r were 
identical, and samples F4, E5, and E3 were also identi- 
cal at all three regions (Figure 5); this finding was not 
surprising given that aza-C r is the result of the wt-FV3 
strain treated with azacytidine and does not necessarily 
represent a strain with a separate evolutionary history 
[22]. It also implies stability in the repeat regions, as 
these regions have not changed between the two strains 
since their initial separation. However, the finding that 
viral isolates from the same geographical area have vari- 
ability may limit the use of STRs as a geographic 
marker. This STR analysis allowed us to better quantify 
the small scale genetic variation that is present in 
highly variable genetic sites amongst FV3 isolates, thus 
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furthering our understanding of genetic variation be- 
yond coding regions. 

The STR analysis we performed in our study has the 
potential to contribute to our understanding of FV3 
tracking and strain designation. Surveillance and phylo- 
geographical analysis of FV3 are pivotal in understand- 
ing how the pathogen varies between different habitat 
sites and amphibian species, as well as for revealing pos- 
sible sources of a disease outbreak [31]. It can also have 
direct effects on conservation by aiding in strategy de- 
velopment to minimize die-offs in high-risk areas, and in 
creating vaccines through knowledge of the FV3 genome 
itself [14,31]. However, exact taxonomic identification of 
viruses in amphibian populations has been difficult given 
the lack of detailed molecular data on FV3 and other 
ranaviruses. The methodology used to classify these vi- 
ruses in the past has been through comparing the major 
capsid protein (MCP) of different viral isolates [3]. How- 
ever, the use of the MCP as a tool to distinguish between 
different ranaviruses, as well as between different strains 
of the same virus, has been called under scrutiny 
[24,32,33]. Thus, there is a need to develop new methods 
of strain tracking for ranaviral isolates. 

The use of STRs for ranavirus strain identification has 
precedence in other virus studies. In one such study, 12 
isolates of human cytomegalovirus (HCMV) were iso- 
lated from various individuals infected with the virus. 
The isolates were then tested for variable repeats in 24 
polymorphic regions, and based on this analysis, each 
viral isolate was designated as an individual strain of 
HCMV [34]. Many of the HCMV repeats used in this 
study were found in non-coding regions of the genome, 
similarly to the ones used in our study. The study sug- 
gested that these changes in repeats are evolutionarily 
neutral and so appropriate for strain identification, not 
only in HCMV, but in other similar, large genome DNA 
viruses [34]. 

Other studies have used coding instead of non-coding 
repeat regions to identify viral strains [35,36]. In our 
study, Region 1 is found in the 19R ORF, unlike the 
non-coding areas of Regions 2 and 3. There are many 
examples of functional microsatellites that are known to 
affect viral characteristics based on copy number, includ- 
ing hepatitis C virus and vesicular stomatitis virus [35,36] . 
Therefore, in addition to being potentially useful in viral 
tracking, STRs from the Region 1 coding region may have 
functional significance in FV3. 

Conclusions 

In this study we have been able to uncover fine scale 
genetic variation between closely related ranaviral iso- 
lates that have different levels of virulence. We have 
shown that substantial genetic variability is present be- 
tween closely related FV3 isolates, both in terms of 



deletions/insertions, and even more so at select STR lo- 
cations. These genomic areas with deletions/insertions 
may present regions that affect viral infectious pheno- 
type, and therefore require investigation. Furthermore, 
we have identified STR regions that may prove useful in 
future phylogeographical tracking and identification of 
ranaviral strains across different environmental regions. 
As FV3 leads to more unexplained lethal infections in 
amphibian populations, studies such as this are neces- 
sary. The genetic insight that they provide into ranaviral 
genomes will prove invaluable when seeking to explain 
the variability of FV3 infectious phenotype in wild popu- 
lations, and in preventing the devastation that ensues. 

Material and methods 

Reagents, viruses 

FV3 isolate SSME, wt-FV3, and aza-C r , were analyzed in 
this study and were provided by Professor Gregory 
Chinchar from the University of Mississipi Medical Center 
(Jackson, MS, USA). SSME was isolated from a wild, spot- 
ted salamander population in Maine, USA, while aza-C r 
was derived from the laboratory wt-FV3 strain through se- 
lection with azacytidine [22]. Both wt-FV3 and aza-C r 
have been previously sequenced and found to be identical, 
with wt-FV3's nucleotide sequence deposited into Gen- 
Bank [12]. Amphibian renal cells (A6 cells) were supplied 
by Niels Bols of the University of Waterloo and main- 
tained in Leibovitzs L-15 media (Invitrogen, Burlington, 
ON) supplemented with 10% fetal bovine serum (FBS; 
Invitrogen, Burlington, ON), penicillin (100 U/ml), and 
streptomycin (100 (ig/ml). 

Infection of tadpoles with FV3 

In order to monitor survival of FV3 infected tadpoles, 
we obtained Lithiobates pipiens tadpoles, approximately 
Gosner stage 25 [37], from the Environment Canada 
Atlantic Laboratory for Environmental Testing in Moncton, 
NB., courtesy of Paula Jackman. The animals were re- 
ceived two weeks prior the beginning of the experiment 
and were kept in 20 L tanks filled with 10 L of aged 
clean dechlorinated water. Each treatment group was done 
as 2 replicates of 25 tadpoles per treatment. Groups tadpoles 
were placed in dechlorinated water, with the host density 
(number of tadpoles per volume of water) adjusted to 1 tad- 
pole per 250 mL to avoid any effect of density on tadpole 
development [38]. To infect, 25 tadpoles were placed in 
50 mL of infected water containing 10,000 pfu/mL of a 
FV3 strain (SSME, wt-FV3, and aza-C r ). According to 
past experiments, such concentration is known to in- 
duce sublethal effects in these laboratory conditions 
[38]. Control individuals were placed within 50 ml of 
FV3-free water. The tadpoles were left within the in- 
fected solution overnight (12 hours) tadpoles were then 
transferred together with the contaminated water in 2 L 
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plastic containers filled with 1 L of dechlorinated water 
(aged for three days) for the rest of the experiment. 
Containers were held in a climatic chamber (Thermo 
Incubator Model 3740) where the temperature was set to 
remain at 22°C with a 12 h:12 h darlclight cycle. Tadpoles 
were fed on a weekly basis after the water was changed 
with standard tadpole food (Carolina Biological Supply 
Company, Burlington, NC) at 30 mg/tadpole for week 1, 
60 mg/tadpole for week 2, and 120 mg/tadpole for week 3 
until the end of the experiment [38]. Starting on week 3 
the water in each tank was replaced once a week with 
clean dechlorinated aged (24 h) water. As a result, exposed 
tadpoles were held in virus-containing water for 3 weeks, 
a period which is long enough for tadpoles to be in close 
proximity with residual infection [38]. Tanks were moni- 
tored on a daily basis. Dead tadpoles were removed to pre- 
vent any scavenging, and stored at -25°C in individual 
plastic vials with ethanol for subsequent analyses.. The ex- 
periment terminated when all the individuals died or 
reached metamorphosis. The procedures used in this ex- 
periment follow protocol #2010-04-02 approved by the 
Laurentian University Animal Care Committee. 

Screening of tadpoles for FV3 

In order to check for ranavirus infection, all animals 
(including euthanized ones) were dissected to remove 
the liver that was then crushed into a 1.5 ml Eppendorf 
tube. The resulting tissue mixture was used for DNA 
extraction. DNA was extracted using QIAmp DNeasy 
Kit following the standard protocol (Qiagen). After ex- 
traction, a double blind PCR was performed using a 
primer set known to successfully amplify a portion of 
the major capsid protein within the FV3 genome: 
MCP-ranavirus-F (5 -GACTTGGCCACTTATGAC-3 ) 
and MCP-ranavirus-R (5 - GTCTCTGGAGAAGAAGAA), 
following the PCR conditions listed in Mao et al. [33], 
using 1.5 ul of template DNA and cycled 40 times. Indi- 
viduals showing two positive amplifications for both 
PCRs were considered infected. We analyzed host sur- 
vival using a survival analysis and failure time analysis 
following the Kaplan & Meier product limit method as- 
sociated with Chi square and Gehans Wilcoxon tests 
(multiple and two sample comparisons respectively) 
[39]. Individuals surviving to the end of the experiment 
were censored to account for our lack of information 
about their true time of death [40] . 

Viral DNA isolation 

The SSME viral isolate was propagated on a confluent 
monolayer of A6 cells, and grown in a 75 cm 2 flask, 
with cells infected at a multiplicity of infection (MOI) 
of 1 PFU/cell. The cells were harvested 5 days post in- 
fection (once cytopathic effect appeared), and viral 
DNA was extracted using the Purelink Viral RNA/DNA 



Mini Kit according to manufacturer s protocol (Invitro- 
gen, Burlington, ON). 

Viral genome sequencing 

Standard kits and protocols developed by the manufacturer 
were used to sequence the SSME sample on a 454 GS-FLX 
platform (Roche Diagnostics Corporation). Briefly, a Rapid 
Library Preparation Kit (Roche, Mississauga, ON) was used 
to mechanically shear 500 ng of template DNA into short 
fragments. A universal sequencing primer that included a 
short DNA sequence unique to the sample (MID tag) was 
then annealed to both ends of each DNA fragment. A GS 
Junior Titanium Emulsion PCR Kit (Roche, Mississauga, 
ON) was used to amplify the sample library, which was 
sequenced using a GS Junior Titanium Sequencing Kit 
(Roche, Mississauga, ON). In order to assemble a full 
genomic sequence, the short sequences produced by 
454 sequencing were aligned with the reference FV3 
genome, wt-FV3, using GS Reference Mapper (Roche, 
Mississauga, ON). Any gaps in the assembled genome 
were then sequenced using custom PCR primers spe- 
cific to each gap, with sequencing performed by the 
Robarts Sequencing Facility (London, ON). The final 
genomic sequence was deposited in GenBank accession 
number KJ175144. 

FV3 sample collection 

FV3 environmental samples were collected from frogs 
caught by hand at various sites along a lakeshore in 
Manitoulin, Kagawong, ON, Canada (Latitude: 45.86418, 
Longitude: -82.27150). The frogs were caught using dis- 
posable gloves which were changed between each animal 
inspection. This method is preferred to the net-catching 
method' as it has been suggested that cross contamin- 
ation can occur via the net. Each individual was toe 
clipped following the protocol #2009-03-04 approved by 
the Laurentian University Animal Care Committee for 
tissue sample collection. DNA was then isolated from toe 
clippings using the DNeasy Blood & Tissue Kit according 
to the manufacturers protocol for total DNA extraction 
from animal tissues (Qiagen, Mississauga, ON). Samples 
then underwent PCR with primers designed to amplify 
specific repeat regions identified in Eaton et al. [19]. 
Primers used included: Region 1-F: CGTGGTCAGACT 
GGTCCTCG; Region 1-R: CACCTCTGTCTCTGAATC 
GG; Region 2-F: GAGTTTACTTGGTGGCCATG; Region 
2-R: TCCTGTCAAGAGATCCCCTC; Region 3-F: CTTG 
CTGCTGCCGTTCAGGC; and Region 3-R: AGAGTGA 
AAAAGGTAAAGGC. 

Sequencing repeats and confirming 454 sequence reads 

In a PCR reaction tube the following reactants were 
combined: 10X PCR buffer (Invitrogen, Burlington, ON), 
50 mM MgCl 2 solution (Invitrogen, Burlington, ON), 5X 
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TAQ DNA polymerase (Invitrogen, Burlington, ON), 
10 mM deoxyribonucleotide triphosphates (dNTPs), 

0. 1 mM primer, 2.5 ng DNA, and water to a final volume 
of 50 ul. The reactions were then placed in a thermocy- 
cler under the following conditions: 94°C for 3 minutes, 
94°C for 30 seconds, 56°C for 1 minute, and 72°C for 
1 minute for 30 cycles. Sequences of PCR products were 
determined by Robarts Research Institute DNA Sequen- 
cing Facility in London, ON, and were analyzed using 
BioEdit v7.0.5. 
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